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ABSTRACT 
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Data from Apollo lunar bistatic radar experiments have been processed to 
give probability density functions for surface slopes. These show best 
agreement with a Hagfors scattering law, though data having both gausslan 
and exponential characteristics also exist. Surface roughness estimates 
range from 4* In maria to at least 8“ In highlands, values which are appropriate 
to 25 m horizontal scales and which are areal averages over tens of square kilo- 
meters. Roughness varies with wavelength, most strongly In maria. 
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INTROnilCTlON 

Dual froqucncy blstatic-radar obst'rvat ions of the moon have been made 
using the Apollo lA, lb, and 16 spacecraft in lunar orbit. Transmissions 
f;*om the conmuind-servico modules were received on earth after reflection 
from the lunar surface. Modulation on the signals, which resulted from 
interaction of the electromagnetic wave with the surface, can be analyzed 
to give information about the latter's properties. In this report we discuss 
prob?bi 1 i ty density functions for surface slope which have been obtained 
from these data. 

The geometry and operational aspects of this expei iment have been dis- 
cussed in detail elsewhere (Tyler and Howard, 1973). Theoretical background 
on the radiowave scattering problem lias l»een developed by Tyler and Ingalls 
(1971). Methods for obtaining the slope density functions from the resultant 
radar echo spectra have been outlined by Parker and Tyler (1973); Lipa and 
Tyler (1976) have applied these technigues to Mars radar data. 

Apollos 14 and 15 provided data at both 13 and 116 cm wavelengths; only 
13 cm was available from Apollo 16. Slope density functions were obtained 
on two degree longi tudinal centers along approximately eguatorial ground bracks. 
The experimentally derived functions have been compared with theoretical 
expressions and have been examined for wavelength dependence. A brief review 
of the background and a summary of the results follows. 
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EXPERIMENT 


In a bistatlc-radar experiment, such as the one conducted here, the 
signal originates at a spacecraft. moving with velocity The incoming wave 
interacts with the surface, which behaves like a polished (but not perfectly 
smooth) sphere. Scattering at radio frequencies is locally specular so that 
most energy reaching the receiver arrives from a region near the point on the 
mean sphere where angle of incidence equals angle of reflection. As the 
spacecraft moves, so does the primary reflection region and over a period 
of time a track is defined across the surface. Spacecraft motion also causes 
the echo to be broadened because of differential Doppler offsets at each 
reflection point within the primary region. 

Our scattering model is based on the assumption that the echo arises 
from specular scattering by those parts of a gently undulating surface 
which are oriented so as to reflect energy toward the receiver. The dimensions 
of the scattering elements are believed to be in the range of one to 1000 times 
the radar wavelength (x), with an effective size being about 200x (Tyler et 
al . . 1971). The size of a Fresnel zone (the area over which a reflection ma ' 

oe considered coherent) marks the upper limit of this range; at the small 
end, scattering passes into the Rayleigh regime. The total echo is the 
vector sum (to account for both amplitude and phase) of the individual 

contributions. Its frequency dispersion is a measure of the surface roughness. 

Scattering facets away from the specular point will contribute less in 

proportion to their probability of being oriented for mirror-like reflection. 
Each contribution to the echo will be Doppler shifted in frequency by an 
amount which is proportional to the scalar product between the look vector towa 
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the facet, u . ani V . Loci of constant Doppler shift arc found to bo along 

d SC 

the intersection of a cone defined by u ,V ■ constant and the planetary 

a sc 

sphere (Tyler and Ingalls. 1971). One can estimate the received signal at 
any frequency by combining this geometrical knowledge with a postulated 
distribution of reflecting facets. 

For this work we will describe facets in terms of their tilt 0 with 
respect to the mean spherical surface. The angle o is indeperdent of azimuth 
and is give by, o ■ arccos (n-ii^) . where h is the unit normal for a facet and 

u^ is the unit normal (radial) to the mean sphere. Any facet for which the 
above condition holds becomes part of the distribution p(0) but only those 
facets for which angle of incidence equals angle of reflection contribute to 
the echo. 

The expression for power spectral density in the received echo w(f) 
may be separated into two parts if one makes reasonable assumptions (Parker 
and Tyler, 1973). The function H(r,s) depends solely on experimental geo- 
metry, while p(0) describes the distribution of roughness. Together, they 
give 


w(f) - yHCr.s) p(0) do 

where the integration is performed over that part of the surface which is mutual Ij 
visible from transmitter and receiver. The vector F is the position of the 
spacecraft and the vector ? is the position of the specular point. 

If a data spectrum W(f) is given, it becomes possible (after invoking 
appropriate assumptions) to invert the computation to obtain P(0), an experi- 
mental slope distribution. In performing the inversion, one assumes that 
the scattering surface is homogeneous over the strips which contribute to 
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Individual frequency bins; there are no a priori assumptions about the 
form of P(0), however, so the result should be an accurate description of 
the surface. In our analysis we have compared experimentally determined 
P(0) with theoretical expressions p(o). 

ANALYSIS 

The inversion techniques of Parker and Tyler (1973) were 
applied to 13 cm wavelength data from each of the three Apollo spacecraft 
and to 116 cm data from Apollos 14 and 15. The resultant slope density 
functions have been conipared with three theoretical expressions. Background 
and Information on their applicability are contained in the Appendix. 

The Gaussian law 

Pq(0) " (sec^O/tan^O^) exp(-tan^O/tan^O^) 

2 

is mathematically convenient, allows for calculation of tan 0 in closed form, 
and has been widely used in theoretical treatments of scattering (see 
Beckmann and Spizzichino 1963). 

The exponential law, 

p^(o) * A exp(-tanO/tanOp) 

does not have the mathematical advantages of the gaussian, but its intrinsic 
properties make it a good alternative to the first distribution. Its use 
in scattering studies follows from substitution into the results of Kodis (1966) 
The third law, derived by Hagfors (1964), 

p^^(O) = B COS0 (cos^O + C sin^O)'^^^ 
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h.is Ihi'ort'Kcal limitations (Rarricl, 1970) but ayroos woll with pla>u»ta»'> 
radar oxpi*rinK*nts (Ivans and Hagfors, l^oR). Points contral to its deviation, 
a gaussian surface height distribution and an exponential surface anticor»*ela- 
tion function, appear to hold for at least some earth surfaces (Hayre and Mi>ore, 


1 % 1 ). 


Our comparison of data with the- above theoretical 
throuvjh consideration of squared error 


expressions was simply 
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K' is computed analytically for best fit, and one adjusts the characteristic 

2 

surface »*oughness (tam'^ or C, depending on the law) to obtain minimum e . 

RESULTS 

Results of the comparisons are shown in Fig. 1 for Apollo 15 data at 

2 

13 cm wavelength. Fig. la shows ) for the three theoretical density 

functions while the lower plot gives the mis surface roughness which one would 
derive from the best fit Hagfoss distribution p^^ (O) (see Appendix). For the 
mi»st part, the Magfors law (solid line in I ig. la) provides the best agreement 
with the data, beyond about 50“ longitude (both east and west) all three laws 
show noticeably higher fit error because of noisier data. 

In maria the Hagfors and exponential laws, both of which have non-negl igible 
tails, are clearly superior to the gaussian. In highlands, the situation is less 
clear, with all three laws showing similar fit errors. From a subjective point 
of view, the fits in highlands do not appear so good as those in maria; the 
normal i/ation procedure used (see Appendix) causes the numerical values 
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for square error to be virtually the same, however. 

The upper limit of roughness values (about 8°) shown for highlands in Fig 
1b may be an artifact of arguit»ent truncation in the experimental slope density 
functions; only probabilities for slopes between 0 ■ 0’ and 0 *20“ were 
obtained in the inversion process. It may be that highland areas have rms 
slopes greater than 8", but that these cannot be determined unless the range 
of 0 is extended. Certainly significant tilts at angles larger than 0 ■ 20“ 
would materially alter the rms surface slope estimate. 

The apparent indifference of the fitting procedure to scattering law in 

highland areas may also be an artifact of the 0 truncation. For the larger 

valu’S of rms roughness, the tails of all dcnsitv functions fall outside the 

0 ■ 20“ limit. Since the absence of a tail appears to be what distinguishes 

the gaussian from the Hagfors and exponential laws in maria, the more favorab1( 
2 

values of c for the gaussian determined using the broad density functions is 
not surprising. 

Fig. 2 shows examples of data which match, under the truncation condition; 
noted above, each of the three theoretical expressions considered. Error bars 
refer to uncertainty in the inversion process; solid lines are best-fit theore- 
tical curves. For the most part, the Hagfors law was the preferred theoretical 
expression, but, as can be seen, the exponential (and occasionally the gaussiar 
sometimes provided the best fit. 

The data were grouped into the broad classifications of maria and highlanc 
for further study. This differentiation was based on a subjective study of 
very low resolution lunar photography; mountainous terrain was considered 
highlands and the large plains were labeled maria. Since each data point 
represents an average response over two degrees of longitude along the 
groundtrack, it was impossible in some cases to arrive at an unambiguous 
classification. Scattering areas which included both mountains and plains or 
which could not bo simply classified as one or vhe other, were labeled 
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“transition'' regions. The breakdown for data obtained along the Apollo 
15 ground track is given in Table I. 

Data points were then plotted on a three variable “phase diagram" 
such as that shown in Fig. 3. The distance from any vertex depends on 
the error incurred during fitting of the data points to the theoretical 
curve. For small Hagfors errors (E^^) relative to exponential and 
gaussian, for example, will be less than either Dp and Dq and so 
the data point will be closest to the Hagfors vertex of the triangle. 

In this way, one can determine the relative importance of Hagfors, 
exponential, and gaussian characteristics in any given set of data. 

The actual Apollo 15 13 cm points are given in Fig. 4. One 
should note the high percentage of points within the Hagfors third 
of the phase diagram. Also, in two-way comparisons, such as the histo- 
grams along each of the edges illustrate, the Hagfors law shows a 
definite preference over either the exponential or the gaussian. Along 
the bottom edge, where the gaussian and exponential lav/s are compared 
independently of the Hagfors, the exponential is favored. 

By separating the data into terrain types, as per Table I, one 
can generate phase diagrams and histograms for mare and highland points 
separately. These show an extreme prejudice against the gaussian law 

when the mare data are used exclusively (Fig. 5). In fact when Apollo 16 mare 
points are plotted in the same way (not shown here), no data appear on the 
left hand side of the diagram at all. The distinction betv/een Hagfors and 
exponential properties in maria is less clear, though there is some skewing 
toward the former. 

When highland points are used exclusively, the preference is toward 
Hagfors characteristics in two-way comparisons (Fig. 6), but the trends 

ORIGINAL PAGE IS 
OF POOR QUALffY: 
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are less strong than with mare data. On the gausslan-exponentlal side of 
the diagram, there Is no clear pattern. There are fewer highland points 
than mare points, so the significance of these latter conclusions must be 
less. Further, the scatter of the highland points suggests that many more 
would be needed In order to draw strong conclusions. The entire highland 
analysis Is Intertwined with the trun.^clon problem discussed earlier. It 
may be that with better and more abundant data the apparent drift from 
Hagfors toward gausslan properties In highlands will turn out to be simply 
an artifact of the 20® cutoff In probability density function determination. 
More data are clearly needed before these questions can be finally answered. 

One should note In Figs. 4 and 5 that there are very few data points 
within the Interior triangle defined by the gausslan and exponential vertices 
and the midpoint of the triangle. This suggests that the Hagfors law Is 
Indeed a good compromise between the gausslan and exponential versions. The 
rounded peak and the non-negl Igible tails which characterize the Hagfors law 
make It good for transition between the two more extreme expressions. The 
fact that most of the data fall within the "compromise" range rather than 
toward either or both extremes Indicates that we have adequately covered the 
range of scattering law possibilities. 

Though not shown here, the Apollo 14 and 16 13 cm data (of lower and 
higher signal-to-nolse ratio, respectively, than Apollo 15) are In excellent 
agreement with these conclusions. 

The VHP data have been analyzed In exactly the same way as the 13 cm 
results. Fig. 7 shows the logarithm of squared error and estimated rms 
surface slope for the Apollo 15 ground track at 116 cm wavelength. Rms 
surface slopes In the v/estern maria are typically 2®, while those In the 
highlands range to 8®. It is this same 8® apparent upper limit in highlands 
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which caused concern regarding a possible truncation error In the processing 
of the 13 cm wavelength data. We thus view conclusions based on these high- 
land results with less certainty than those for maria. Tyler (1978), on the 
other hand, has fourid both theoretical and experimental evidence for believing 
that the change In be/:av1or between maria and highlands Is real. 

The surface roughness estimates from the 13 and 116 cm data sets (Figs. 1 
and 7, respectively) have been compared In Fig. 8. EarMer studies (Hagfors 
and Evans, 1968; Parker and Tyler, 1973) have shown wavelength dependence In back- 
scatter and, from a ‘limited analysis, wavelength dependence In oblique scatter. 

The Inference from Fig. 8 is that there Is a strong wavelength dependence In 
the scattering process thoughout mare regions. In the western maria the factor 
of two difference In roughness between 13 cm and 116 cm wavelengths Is consis- 
tent with a surface roughness (s' variability of form 

s « X 

where the exponent a has a value on the order of -0.3. In Montes Apenninus 
there Is no apparent wavelength variation; roughness at both wavelengths 
appears to be about 6-8* so that a In the above expression would be 0.0. Mare 
Serenitatls, which appears to have wavelength dependence, but to a lesser 
degree than the western maria, shows a dependence of approximately a ■ -0.2. 

A tabulation of the wavelegth results by geographic area Is given In Table II. 
One-sigma error bars represent the variability within the N data points and 
not the uncertainty In obtaining those points. 
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SLOI»K DISlHlIUrriON FUNCTIONS ANH R/\DAH CROSS SKCTION 


Several commonly used probability density functions for surface slopes 
are given. Both unidirectional slope (slo|)e measured along a particular 
line, such as the x-axls) and adirectlonal slope (the title of a facet with 
respect to the normal to the mean surface) are used. 

1. ReOj Hontherington , and Mifflin have defined tilt through 

where dA is the area of an incremental scattering facet tilted an angle 9 
with respect to the mean surface normal, and du) is the element of solid angle 
into which the normal to dA points. The quantity is the density func- 

tlon giving the probability that sold facet normal is within dco. 


Projected area on the mean surface is 


dS B CO80 dA B COS0 Po,,,,(0) 

HHM 


and the normallxatlon requirement is 


/dS . f.. 


fcos0 dA B 

hemisphere hemisphere hemisphere 




Ah thv clunnlty function 1 h 1 ndopi’ndi'nt of nzlmuthnl anKlc, Uio no rinn 11 motion 
may be rewritten 


n/2 

f d0 ■ l 

0 

where du * sing d9 d0, 

Rea, Hentherington, and Mifflin ( 196 ^) modified their work somewhat 

after ItagforH pointed out that Pu,„,(6) was defined with respect to 

KHM 

the actual undulating surface rather than to the meon flat surface, A new 
density function can then be defined 


p„(e) . co«e 

with normalization 


.p / 


Pjj(e) sine dO 


2. Muhlema n ( I 96 IO defined the probability density function Pj^jCSj®) 
whore Pjj(0»0) dio is the probability that a facet normal points Into solid 
angle d'J) at zenith-azimuthal coordinates (0,0). The proper normalization 
in this case is 


2n Jt/2 

If 


0 0 


Pm( 6,0) sln0 d0 d0 


1 


ORIGINAL PAGE IS 
OF POOR QUALITYi 


ttttually litdi'ptMuli'nct* of (1 Ik iiKKUir.cil ku that 


r ./2 

^ Bine c10 . 1 

0 

IMrki’r (lf’^ 3 )» foHowlnc MuMoman, introihictU a onr-dimi'niiloitnl 
dfUBlty function 


Pp (e) 


/ 


p „ ( e , o ) do 

m 


where norma llr.nt Ion 1 h nccordlim to 


/ 


p^^(e) sine de 


1 


3. Beckmann (l9t>3) » nmonp others, hog an 
for the rouKh surface J(x,y) from which a dcnsi 
terms of partial derivatives can be ol)tnincd. 


with an ana ytlc expression 
ty function Pj,(C^.Cy) *n 
Normalization is through 


CD OD 

/ f »1. "tx 

_CD — CU 


e 1 


1» , Slo pe proba bility density has been described by Hagfors ( I 968 ) 
through thf‘ function p (s,0) where s c tan?. This expression is normal- 
i ze<l using 



Pj. ‘'S dO sr 
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whoro convorNlon to oy Itiutrlcnl nnt i-k 1 h imuU- )Htshll>tf thrmiKli 


2 

H 



Tlu’ tU'iiHity it* usually Imlt'poruU-i.t of It iiml out* iucium n I ly pl■^•fl'r» 

to uso 


/ 




*). Ctimparlsou of tho abovo functions can boj;ln by equating; the* 
integrands of S€>ctlons 2 and l| above 

Pnj(0,O) slnO d0 dO a p^(s,0)s ds dO 

m p (tan0,(l) lanO sec 0 d0 dO 

This leads to 

■ sec e p (tanO, O) 


From sections P and It wo also have 


'V 
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Pj^,(e,o) do) 


0 


OUK'.VNM^PA^ 

OV' l>00ll Ul3AUt» 


P„, (tane) 


f. 71 

/ 


P^ (tanO.tt) dtl 


and 


HO that 


P„ (0) 


Pn P„(0.t’) 


m I'n 0 p^(tnn0,(O 

H H»'i‘ 0 p , ( t nnO ) 


From HoctloHM 1 nml P wt* have, roRpoot 1 vrly , 


n/P 




) hIiiO (Id m 1 


0 


and 


n/2 

/ 


P^(0) Hind (Id 


. 1 


From these om* oMalns 


Pp(0) - •’« P„(S) 


and thus 


P|,(0) ■ 0 P^(t:>»0,l’) 


sec'^0 


n K 


p /tnnO) 
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Finally, from ni'otlitn 1, »«* luivi* 


P„(o) - co.e i<s) 




;’n " ‘V'(tiinO) 


Wo can nlHO equate the Intei'ramlx of HOcltonH 3 ”'*‘l *•» 


P (C ,C ^ ‘>C ‘>C 

'H''’x*'’y ''y 


p (s,0)s ils ilO 


p (tnnOjl') tanO Kec 0 ‘10 ‘h’ 


In this case, the conversion between ,C ) coonllnales ami (s,0) 

X y 

coordinates is r t ralpht forward and 


where 


p^(tan0,tl) « l’n(w^«Cy) 


Cx' - S‘ 


tan 0 


6. Hatinr croHn m’o t Io n can be* obtained In nrvcral wuyii. Fnl lowing 
the notation of hc*i*lon I, we can dtnote |K»wer Intercept «>d by a facet n» 


dA COM 1 


wheie Is power density In the Incomlni; plane wave and 1 Is nii|;le of 
Incidence on the facet. 


This pt>wer In reflected Into the forward direction with density 
P^. lUit direction of the normal to dA is uncertain by an amount dn so the 
direction of the reflected ray(s) Is uncertain by 


du m It cos 1 d'u 
r 


At a distance r from the facet, that means ix»wer can be spread over an 
sr ’ 


di' r 
r sr 


The expected power density at a receiving antenna, If n Is uniformly 
distributed over d» and Kcometry is not otherwise too perverse, Is 


Pj dA cos 1 

dn r ^ 
r sr 


The radar equation Rives expected power from a recelvlnn antenna of 
aperture A^^ as 


hnr 

sr 
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I’^ ilA con i 

di> r *’ 
r nr 
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|{ 
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•ifl I'Oh 


dA 
d n 


Substitution for dA anti di'^ bIvoh 

"o ■ ” "kihi («’ 

whore thi* convornion to ra«lar cross section per unit area of the mean surface 
is strnipht forward. 


Kadar cross section 

1* iMI 

also be expresseil in terms c»f other density 

functions as 



a 

o 

n 

It 

n sec 9 p^ (tan 0,0) 

0 

o 

o 

It, 

P Pj,' (‘«'» 0) 

tr 

o 

C7 

n sec 0 P|| (e) 

a 

o 

xs 

® P„ (el 

c 

o 

ts 

" ® ”11 

The last expression above 

' Is 

identical to th«> one derived by Ilarrick (1969) 


for scatterlnR from a gently umlulatini; surface. 


7. Ci aiiHHlun Surfacc M nrc among Uio rnoHl cwiHlly hntuilpd mathematically. 
Deckmann (I 963 ) u»e» a gnuHRlnn height dlstrlhutlon ond gminslan surface 
auto-corrclntlon function to arrive ut a slope density function of the form 


o o 


where Q Is the unidirectional rms slope. This converts directly to 
o 


(tan 0, 0) 


exp 


2 

tan 0 


2«C, 




and thence to 


Pg, (t«'» 0) 


I ( tan^0 1 

r ^c. = I 


From section 5 these can be 

( 0 ) 

P,l (0) “ 

Pp (0) - 


rewritten to 


sec 9 


exp 


3 

secj0 

2nC 
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exp 


3 

sec 0 


exp 


give 
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tan^0 


I _ tnn^9 
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tan 9 
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i ec. 
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I 
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Kadnr croRH Rcctlonii per unit area for the above density functions are 


KHM 


U (?) 

sec 6 I tan 6 ( 

2 ■ 2 


- ( 9 »®) - ( 9 . 0 ) - O (e) - a (e) 


\ *ec 6 

- o ®’«P 

2 C_ ^ 


Op'^x' y' 


'll 

l.iilliZ 
I 2 C ^ 

X„ 


Dual argument cross sections ore those derived from density functions 

requiring normalization over two variables; single argument cross sectioi s 

are those from density functions requiring normalization over only one 

variable. In either case, cr must be integrated over all space if total 

o 

scattered power is desired. 


9. Hoot mean square slope can be computed in a number of ways. The 
easiest Is by using 
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t1/2 


(s) s ds 


For the density function p^, (s) given in section 7, this lends to nn 
ndircctionnl rms slope, or tilt, of 


s 

o 
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Kor P.irkff** dvniilty function, the ndl rvet lonu 1 rra« k1o|h> Is obtained 
UNlnR 


s 

u 



tail 0 Pp(®) "InS ‘IS 


1/ii 


One also obtains 


s 

o 



when the function of section 7 Inserted, 


For the other density functions, ndlrectlonal rms slope can be 
found from 


r 2n 


T 1/.? 


// 


o o 




it/2 1 1/2 
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2n / tnn 0 


/ ■■ 


n/2 
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tan 0 Pjj(0) 'ln0 dO 
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1/2 


r 2n n/2 

// ■ 

L O O 


an 6 Pj^j(^»^) d0 dO 


1/2 
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o 




♦ c,'’) "„(c,-Cy) 
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For an Inotropic aurfacc, havlni; the sniiic atntlatica In x and y direction* 

and independent alopes ; and C , above exprenaion* may bo 

X y 

written 


a 

o 



c/ p„, (c,) i;. 


1/2 


where p„. (C ) is the one dimcnaionnl analog of Pn(v„»C„): 
B’ ' x' u X y 


p„(Cx,C ) - 


P|ji(Cx) ^^ y ^ 


In practice, val'ics of radar cross section as a function of incidence 
angle .iro used to determine the corresponding probability density function 
according to che relationship of section 6. Rmn slope then is obtained ns 


in this section above. 
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TERRAIN CLASSIFICATION ALONG APOLLO 15 GROUND TRACK 
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These data groupings are subsets of the N = 27 Mare Imbrium and Mare 
Procellarum group immediately above 


FIGURE CAPTIONS 


1. (a) Fit error between experimental surface slope probability 
density functions and gaussian (+), exponential (x), and Hagfors (solid line) 
theoretical curves. Hagfors is somewhat better than the exponential; both 
are considerably better than the gaussian. (b) RMS surface roughness approp- 
riate to best fit Hagfors law along Apollo 15 ground track at 13 cm wavelength. 

2. Examples of good fits along the Apollo 16 ground track to each of 
the three theoretical expressions considered. Error bars refer to uncertainty 
in the data points and not to the fit between data and theoretical curves. 

3. "Phase diagram" used to show data behavior with regard to each 

theoretical scattering law. Distance from a vertex is based on the sqjare 

error in the fitting process; the constant K is chosen to satisfy the identity. 
2 2 2 

In the figure, is less than either Cq or (indicating that the 

Hagfors law provides the best fit) so that D„ is less than either or 

n b 

and the point is closest to the Hagfors vertex. 

4. "Phase diagram" showing all 13 cm wavelength Apollo 15 points. Most 
surface slope probability density functions show best agreement with the 
Hagfors law; a lesser number show best agreement with the exponential, and 
the least number show best agreement with the gau'^sian. Histograms show 
distribution of points when only tv/o scattering laws are compared. 

5. "Phase diagram" for 13 cm Apollo 15 points in maria. When highland 
points are omitted, the distribution moves away from the gaussian scattering 
law. 



6. "Phase diagram" for 13 cm Apollo 15 points in highlands. Consider- 
ation of highland data only shows Hagfors preferred but with a shift away 
from the exponential and towa d the gaussian. 

7. Fit error (a) and Hagfors law rms surface slope (b) from Apollo 15 
data at 116 cm wavelength. Fit error shows same behavior as at 13 cm (see 
Fig. 1); surface roughness is noticeably lower in maria but similar in high- 
lands. Sharp spike near 30“W longitude results from very irregular terrain 
around the crater Euler. 

8. Ratio of roughness estimates at 116 cm to those at 13 cm along the 
Apollo 15 ground track (see Table II). 
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